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Abstract 

The evolution of disks around Young Stellar Objects (YSOs) is 
deeply affected by the YSOs ultraviolet (UV) radiation field especially 
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in the 500-1100 A spectral range. This two dominant processes are; the 
photo-dissociation of H2 molecules in the Werner and Lyman bands, 
and the emission of photo-electrons from dust grains when high energy 
photons are absorbed. Photo-electrons are an important source of gas 
heating. In this letter, dust grain charging when exposed to various 
possible UV fields in the YSOs environment is investigated. Numerical 
simulation of the evolution of photo-electrons in the electric field cre- 
ated by the charged dust grains are carried out to obtain the charging 
profile of dust grains. Prom the simulations it appears that the differ- 
ent spectra produce significant quantitative and qualitative different 
charging processes. Both the UV background and the Ae-Herbig star 
radiation field produce a relatively slow charging of dust grains due 
to the low fraction of sufficiently energetic photons. The radiation 
field of T Tauri Stars (TTSs) is harder due to the release of magnetic 
energy in the dense magnetospheric environment. These numerical 
results have been used to propose a new simple analytical model for 
grain charging in the atmosphere of protostellar disks around TTSs 
susceptible to be used in any disk modeling. It has been found that 
the yield decreases exponentially with the dust charge and that two 
populations of photoelectrons are produced: a low energy population 
with mean kinetic energy E = 2.5 eV and a high energy population 
with = 5.5 — 6 eV; the energy dispersion within the populations is 
~ 1.3 eV (T ~ 1.5 X 10^ K). The high energy population is suscepti- 
ble of dissociating the H2 and ionizing some low ionization potential 
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species, such as the Mg. These results add an additional role to dust 
on the chemistry of the layers just below the H2 photoionization front. 
This photoelectic yield has been applied to a simple evaluation of the 
dust charge in the atmospheres of accretion disks (a-disks). 

1 Introduction 

UV radiation strongly affects matter and the circumstellar environments 
around Young Stellar Objects (YSOs) are very sensitive to it. Of special 
interest is the interaction between the stellar UV radiation field and the 
accretion disk that channels mass infall onto the star. 

In the diffuse atmosphere of the disk, the dominant interaction is the 
absorption of UV photons by the H2 molecule in the Lyman and Werner 
bands; about 15% of the photons cause photo-dissotiation of the rather inert 
H2 molecule [5] leading to both the formation of a photo-ionization front 
that propagates inwards the disk and an outward photo-evaporative flow on 
the surface, which carves the disk and produces a slow thermal wind [I2] 
that has a fundamental role in disk evolution [1]. Indirect evidences of the 
absorption of UV radiation in young planetary disks come from the detection 
of H2 electronic transitions that are pumped by the Lyman a radiation from 
the star |TT]. 

Though II2 by far is the main component of molecular clouds and proto- 
stellar disks — the gas to dust mass ratio is ~ 100 and CO/H2 ~ 7 x 10^^ — 
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dust plays a relevant role as source of extinction and free electrons, not to 
mention its role in disk chemistry. The largest cross-section of a dust grain 



both towards X-rays and optical wavelengths) thus, dust is more opaque to 
the optical/ultraviolet stellar spectrum than to the infrared radiation pro- 
duced by itself. Furthermore does the interaction between dust grains and 
the UV radiation field result in the generation of a population of high en- 
ergy photo-electrons. There is evidence of such hot electrons co-existing with 
molecular hydrogen around T Tauri stars (TTSs) [H] and in the shocks of 
YSOs outflows with the environment, like in the Herbig Haro Object (HHO) 
HH 2. As an example, the observed H2 emission spectrum of HH 2H can be 
fitted by radiation from H2 at 1000 K that is coUisionally excited by 20 eV 
electrons [18]. There are, however, two problems with this interpretation. 
Firstly, it is unclear how hydrogen molecules can mix with hot electrons and 
radiate without dissociating. Secondly, the H2 temperature must be kept 
above ~ 1000 K to excite the vibrational levels and below ~ 10, 000 K to 
avoid dissociation. Photo-electrons from UV irradiated dust could provide 
the needed thermostat, given that they are sufficiently energetic. Heating of 
the gas by collisions with photo-electrons depends on their emission rate. Jet 
and their kinetic energy through: 



arises when interacting with radiation at wavelenght; 



s ~ 750 A (declining 




(1) 
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where {E — eU) is the kinetic energy of an escaping photo-electron, E is 
the energy of the impinging photons and U the barrier to escape the grain 
potential, uh and a are the number density of H2 and the collision cross 
section of H2 molecules with the photo-electrons, in the simplest version of 
a dust grains fluid embedded in molecular gas. It should be noted that the 
kinetic energy of the photo-electrons depends on their interaction with the 
electric field generated by the charged electric grains and it is directly related 
with the grains charging profile. 

In this paper, the results from detailed numerical simulations of the charg- 
ing process of dust grains are reported, for various UV radiation fields of in- 
terest for star formation research (photoelectric effects by X-rays and Auger 
effects are beyond the scope of this paper) Also, a simple parametrization of 
Jg, the photoelectric current produced by a dust grain, is derived for grains 
submitted to the irradiation of the magnetosphere of a TTS. The dependence 
of the photolectric yield on the dust grain charge is analyzed in Sect. 4 and 
the population of ejected photo-electrons is characterized in terms of it ki- 
netic energy. Finally, the relevance of photoelectric charging of dust in the 
atmospheres of a-disks is discussed. 

2 Numerical simulations 

The applied scheme to account for the interaction between incoming photons 
and dust grains follows four steps: 
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1. Capture of ionizing photons by the dust grain. 

2. Yield of photo-electrons released from the surface of the grain-core. 

3. Determination of the initial kinetic energy of a released photo-electron. 

4. Dynamics of the electrons in the local electric field that results from 
the charged grain-core and the surrounding electronic cloud. 

The first three phases are based on the scheme originally proposed by P] for 
the treatment of the interaction between radiation and dust grains. All inter- 
actions are parametrized for the most common silicate dust grains (see [T71 

m)- 

2.1 Photon capture 

For a photon with the wavelength A, the cross section for capture/absorption 
is given by: 

(^{r grain/ A, '^{u)) = Qabs{r grain/ A, )7rrJ^„i„ (2) 

where n is the complex refractive index for the grain material and r grain /^ is 
the ratio between the grain size and the wavelength of the impinging photon. 
The absorption efficiency Qabs is obtained by applying the MIE-theory using 
the functional form of the refraction index for silicate as stated by [1]. 
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2.2 Yield of photo-electrons 

For a bulk sample the electron yield produced by a photon, Y, is approxi- 
mated following the scheme by [28j : 




where Q = hu — Xmou and the parameters Xmoi = 8 eV , a = 0.5 and /3 = 5 
represent silicates. For this yield model to be valid the impinging photon 
must interact with a molecular orbital and its energy must be < Xatom- For 
the more energetic photons a decrease in yield occurs since they interact with 
the atomic orbitals, which results in a weaker coupling. To account for this 
behavior we suggest a simple and analytic approximation: 



where T = hu — Xatom whereas a and b are fitting parameters. The values 
used for silicates are: Xatom = a = 20 eV and b = —2. For these values a 
good fit to the more accurate model by [21] is obtained, see Fig. [TJ With 
this additional term the full expression for the yield of photo-electrons from 
a bulk sample becomes: 




F(0) 



a(0/XmoO 



(3) 



l + f3{Q/Xmol) 




(4) 
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Size effects result in a dependency between yield and the curvature of the 
grain surface as an increased curvature results in a stronger electric field 
within the grain-core. Moreover, the distance a photo-electron has to travel 
within the grain — to reach the surface — decreases. These effects have been 
modeled accurately by ^27^ however, the approximate model as suggested 
by [3] has been implemented into our numerical code: 

_ /3' a^-2a + 2-2exp[-a] 
nV»nJ ^ ;32_ 2/3 + 2 -2exp[-/3]' ^ ^ 

where (3 = r grain /I photon is the inverse penetration depth of the photon, a = 
r grain/ Iphoton + T grain / 1 electron relates the penetration depth of the photon, 
^photon, and the length for a photo-electron to travel within the grain-core, 
hiectron- This distaucc is assumed to be hiectrmi = 10 A and represents the 
distance from the molecular orbitals containing loosely bound electrons to 
the surface of the grain. The photon penetration depth is given by Iphoton = 
X[A7T'<s{n)]~^ , where A is the wavelength of the impinging photon and Q{n) 
is the imaginary part of the refractive index of the grain material. The total 
yield is then given by: 

Ytot = Y{hu)Y{rgra^n) (7) 

Photo-electrons with a total yield exceeding unity are always be released. 
When the yield is less than unity a Monte Carlo approach is applied and a 
photo-electron only gets released if the determined yield exceeds a random 
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number drawn from a normalized even distribution. 



2.3 Initial kinetic energy of the released photo-electrons 

The distribution function for the initial kinetic energy of photo-electrons is 
assumed to be parabolically shaped as suggested by [28] and a normalized 
distribution might be given by: 



-6^2 + QEE* 
""high 



f{E) = " ^.r (8) 



where E^^^^^ is the excess energy of the impinging photon when the ionization 
energy has been subtracted i.e., the highest possible kinetic energy a released 
photo-electron might get. 

As the distribution is normalized for the interval [0, E^-^i^] a random num- 
ber, fi, between zero and one is used to set the kinetic energy, -Ekin, for the 
released photo-electron through the expression: 



^kin _6E2 + 6EE* 
^^ = I -3 '-^dE (9) 

^high 



where Root() is the root of the polynomial within the above stated energy 
interval. 
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2.4 Dynamics of electrons in the electric field of a 
charged dust grain 

In the event that a photo-electron is released, a stationary positive charge is 
assigned to the center of the dust grain, which is assumed to be a spherical 
core with the radius R. This radius is furthermore used to determine if an 
orbiting electron recombines with the charged grain-core, which is the case 
if it gets within a distance of R from the center of the grain. 

The electron emission site from the grain surface is selected randomly 
on the hemisphere facing the source. The site is obtained from a three 
component vector drawn from a Gaussian distribution scaled by the core 
radius. Under the assumption that the photon-source is placed at x = oo, it 
is only the half-sphere where x is positive that undergoes emission processes 
and the site vector becomes [Abs(rr), r^, r^]. Besides defining the emission 
site the vector also represents the surface normal — the tangential plane at 
the emission site — used to determine the half-sphere, which the direction of 
the initial velocity vector [vx,Vy,Vz] is restricted to. 

While an electron orbits a charged grain-core it might escape the local 
environment that occurs if the boundary of the simulation cell is reached; a 
cube of side length 100 /im. At this point, is the background electric field 
assumed to become dominant over the local field. To simulate the trajectory 
of the released electrons a Velocity Verlet algorithm (sec Appendix A) and 
an accurate Coulomb potential are applied; the time step is 2.5 x 10~^^ s. 
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The grain-core is kept fixed in the center of the simulation celL 

3 Dust grains interaction with YSOs radia- 
tion fields 

Four irradiation fields have been selected for this study. Two of them are 
black bodies at temperatures of 50,000 K and 10,000 K simulating the photo- 
spheric radiation fields from O and A stars, respectively. A third spectrum is 
the fit to the interstellar radiation field made by [5] based on the early mea- 
surements of the Copernicus and Voyager missions [10] . This soft spectrum 
is used as a standard in chemical modeling of the disk; often the intensity 
of the radiation field is scaled with respect to the mean interstellar radiation 
field by means of the standard G factor (typically G = 10^ at 100 AU i.e. 
[20]). Finally, a hard UV spectrum has been built for this study to model 
the hard energy radiative output of the TTSs. Notice that there is little 
information on the spectral energy distribution of the TTSs in the range 
10 A < A < 950 A. Furthermore, only few TTSs (TW Hya, RU Lup and 
T Tau) have been observed in the 950-1100 A range with the Far Ultravio- 
let Spectroscopic Explorer [301 E] • There are however, strong evidences of 
energetic processes ongoing in the atmospheric/magnetospheric environment 
of the TTSs. On the one hand, the near UV (1150 A < A < 3200 A) indi- 
cates atmospheric structures extending to the inner border of the gas disk 
(see [7] for a review). On the other hand, the observed X-ray luminosity 
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of YSOs is on average about lO^^erg and the energy distribution is well 
fitted by the thermal radiation of optically thin plasmas at temperatures of 
~5-30 MK (see i.e. for comparison the Solar X-ray luminosity is about 
lO^^erg s~^ and the plasma temperatures 1-2 MK. Thus, the propagation 
of the X-ray radiation on the dense (> 10^ cm"'^) and extended (2-4 R*) 
atmospheres/magnetospheres of the TTSs is expected to produce a rich and 
energetic spectrum. The spectrum used in this work has been used to fit the 
observed ClV/Si III] and C 111] /Si 111] lines ratios [8]. 

The spectral energy distributions (in photons s~^cm~^) have been nor- 
malized to the same bolometric flux. The ionizing part of the TTSs model 
spectrum and the black body at 50,000 K are plotted in Fig. [21 also the 
cumulative distributions of ionizing photons are plotted for reference. Notice 
that all the energy distributions are smooth besides the TTS distribution. 

The simulated charging profiles of a 30 nm dust grain are shown in Fig. [3] 
for the four energy distributions: black body at 50,000 K (BB-50kK), black 
body at 10,000 K (BB-lOkK), UV background (UV-BG) and TTS. From the 
simulated results, it appears that the TTS spectrum gives raise to a S-shaped 
profile whereas linear profiles are obtained for the two black body spectra. 
The UV-BG saturates at a charge of approximately 125 e. It also appears 
that the quantitative charging when exposed to the TTS or the BB-50kK is 
alike and that the UV-BG or BB-lOkK gives raise to the similar charging 
rates. The relative slow charging obtained with the UV-BG or BB-lOkK is 
a consequence of a low ratio of photons sufficiently energetic to both ionize 
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an atom and provide the needed kinetic energy for this released electron to 
escape the charged grain-core. Whereas the relative fast charging observed 
for the TTS or BB-50kK is possible because of a larger ratio of energetic 
photons. 

The saturation threshold apparent in the figure for the UV-BG spectrum 
is caused by the high energy cut-off of the energy distribution. The shortest 
wavelength corresponds to a photon energy of 14 eV that gives an excess 
energy of 6 eV, which equals the needed kinetic energy an electron requires 
to escape the electric field created by a grain of charge 125 e and a radius 
of 30 nm that explains the observed plateau level. The linear profiles for the 
BB spectra are consistent with the experimental observation by [22] and is 
a behavior that is accounted for by theoretical models for the charging of a 
capacitor. From the linearity of the charging profiles it can be concluded that 
effects caused by the interactions between released electrons are negligible for 
the simulated photon fluxes. The S-shaped charging proflle caused by TTS 
is a consequence of the line structure of this spectrum. From the cumulative 
distributions it also appears that the ratio of photons at high energies is 
larger in the TTS spectrum than in BB-50kK, see Fig. [2l A consequence of 
this shows in the interval from [1 /is, 10 fis] in Fig. [3] where the TTS almost 
manages to reach the same charge state as obtained by the BB-50kK. If the 
energy scale shown rightmost in Fig. [3] is compared with the one uppermost 
in Fig. [2] it emerges that the energy level at which the curvature changes after 
10^^ s is coinciding with the flrst signiflcant peak in the TTS spectrum. 
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3.1 Dependence on the size of the dust grain 

A series of simulations have been carried out to address the dependency 
between the grain size and the charging profile. Only the TTS spectrum 
is used for this purpose. As expected, the maximum charge a grain may 
achieve increases as the grain size increases. However, the maximal charge 
is determined by the radius of the grain rather than its volume, see inset in 
Fig. [3l As a result, the maximal charge-to-mass ratio increases as the grains 
becomes smaller. This behavior is caused by the electric potential at the 
grain surface since the energy barrier an electron must overcome to escape 
the grain surface is oc Q/rgrain- 

4 Model of the photoelectric emission in the 
TTSs environment 

Two main results can be derived from the simulations: (1) the charging 
profile qualitatively remains unaffected when the grain size varies and (2) 
the maximal charge obtainable by a grain is solely determined by its radius. 
Thus, to reproduce the charging profile of an arbitrary sized grain it is suf- 
ficient: (1) to fit the general functional form of the charging profile for a 
grain and (2) multiply it by a scaling factor, which is a linear function of the 
radius. 

The charging profile of a fiducial 10 nm dust grain has been fit. The grain 
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is irradiated with photons from the TTS spectrum until it reaches a satura- 
tion threshold of 275 e. During the simulation, 1-2 electrons usually orbit the 
grain before they either escape or recombine with the grain core. Electrons 
leaving the computational box are lost to the background. The charging 
profile has been normalized to a photon rate, S'19, of 10^^ s~^cm~^ with en- 
ergies above the dust grain threshold for photo-ionization that corresponds 
to photons with wavelengths A < 1550 A. This photon flux is equivalent to 
what is expected from an extended magnetosphere of radius 3 Rq producing 
a total ionizing luminosity of 0.03 Lq. Notice that for the TTS spectrum, 
this photon rate is 4.8 times higher than the photon rate in the Lyman con- 
tinuum (1100 A< A < 912 A). The charging profile is represented in Fig. 4 
and the polynomial fit to it is given by, 

logio Q{x) = -0.0015x^ - 0.00455x^ - 0.0177x^ + 0.0213x + 2.4326 (11) 
. = log:o(5.<) = log,„ (^^j^-^ J-) e 1-8.4, -1.3] (12) 

The functional dependence for Q{x,r grain) is then described by: 

Si9t < 4 X 10"^ 

4 X lO-'' < ^19^ < 5 X 10-2 (13) 
§19^ > 5 X 10-2 

For this grain size the asymptotic value for the final charge is set to 275 e as 



Q{p^i grain) 



^' grain 

10 nm 



if 

Q{x) if 
275 if 
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only photons with wavelengths less than ~ 258 A will be able to cause any 
further charging that are assumed to be of an insignificant rate, see Fig. |2j 
The charging profile for the proposed analytical expression and the simulated 
data are shown in Fig. |H From the figure it appears that the proposed model 
reproduces the main features of the numerical simulations being an abrupt 
initial charging and a final plateau level. For the remaining intermediate 
interval a good agreement is obtained between the simple analytical model 
and the more accurate numerical results. 

4.1 The dependence of the photoelectric current on 
the grain charge 

To derive a simple parametrization of J° on the grain charge Q, a new series 
of simulations was conducted where the grain cores are kept at a constant 
charge level. In these simulations the kinetic energy of electrons escaping 
the grain is recorded however, rather than loosing these electrons to the 
background, they are reassigned to the grain-core to maintain its charge 
state. To generate the data grains at different charge levels are exposed to 
the TTS spectrum for 1 fis. Exactly, the same series of photons is used 
for each run. The number of ejected photo-electrons is found to decrease 
exponentially as, oc iQ-^-^'^^Q^ where Q is the charge of the dust grain and 
the current of photo-electrons, J|, emitted from the 10 nm grain, might be 
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expressed as, 



r = 6.3 X 102 electrons ^^^^^_o.oi6Q (^4) 

/iS 



Under the assumption that photo-electrons are emitted isotropically from 
the 10 nm dust grain surface, the photoelectric flux is given by. 



^1 = 2 X 1020 electrons ^^^^^_o.o,eQ (^5) 
cm^s 



that can be scaled as, 

= 2 X io2o electrons ^^^ / ^ ' ^q-o.oi6Q (^g) 

cm^s VIO nm/ 

for any arbitrary radius of dust grains, where the dependence of the dust 
grain charge on the radius is inserted from EqJT3l 

4.2 Energy distribution of the photo-electrons 

The energy distribution of the ejected photo-electrons follows a bimodal dis- 
tribution with a narrow peak at about 2.5 eV and a broad shoulder peaking 
at ~ 5.5 — 6 eV. As the charge of the grain grows — which results in an 
increased local electric field — only electrons in the high energy tail of the 
distribution are able to escape and the electrons flux decreases dramatically, 
as shown in Fig. [51 This distribution is caused by the energy distribution 
of the TTSs spectrum irradiating the dust grain; the first peak at 2.5 eV is 
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produced by the Lya continuum while the shoulder at 5.8 eV is caused by 
high energy far UV photons produced by highly ionized species in the 700A 
- 850A wavelength range. A useful parametrization of this behavior is given 

by, 

where E' = E + 0.14(Q — 1) and J|(Q, E) is the number of photo-electrons 
ejected with energy in the range {E — 0.5, E + 0.5) eV, for a grain with 
charge Q and radius 10 nm. The fraction of photo-electrons ejected with 
energies susceptible to ionize the hydrogen (13.6 eV) is 3% for Q = 1. The 
parametrization of the photo-electrons energy distribution allows to address 
the effect of the UV irradiation, in terms of the two electron populations 
produced: 

• A low energy population with E = 2.5 eV and dispersion 1.2 eV. 

• A high energy population with E = 5.8 eV and dispersion 1.4 eV. 

where the dispersion is a measure of the temperature of the electron beam: 
1.4 X lO'* K and 1.6 x lO'' K for the low and high energy populations respec- 
tively. Notice that the high energy population is well above the dissociation 
threshold of the H2 molecule (4.52 eV) and may contribute to the ionization 
of Mg {xi = 7.64 eV), that is a fundamental charge carrier for protostellar 
disks [T61II3]. 
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5 Discussion: dust charging in TTSs' disks 

The dust charging profile derived in Sect. 4, provides a useful parametrization 
of the photoelectric current generated by dust grains irradiated by a TTS 
spectrum. To determine the final charge of dust grains is a much more 
complex problem since other processes must be taken into account, including 
collisions with charged particles or charge transfer reactions between the dust 
grains and the surrounding molecules (see [SB] [20] for a detailed evaluation of 
the ionization degree and the main chemical processes in protostellar disks). 
For a generic electron density, ne, and temperature, T^, the final dust grain 
equilibrium charge is determined by the balance of the photoelectric emission 
current, $g, and the driven into the positively charged grain by the collisions 
with the surrounding thermal electrons. The fiux of absorbed electrons per 
unit surface on the dust grain, [23j is given by, 

y-nniej \ TgrainkT J 

where is the sticking probability of the electrons striking the grain and 

/ \ 1/2 

I ) is the thermal velocity of the electrons. For a 10 nm radius dust 



gram, 



e^O 10^ K 

^ i.erg^V^ (19) 



grain ^ 
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and thus, 

Both, $g and $g depend on the environmental conditions. Let us make 
an estimate of the relevance of photoelectric charging in the atmospheres of 
TTSs' disks. 

5.1 Physical properties of a-disks 

Assuming a simple a-disk model [21] , the variation of the disk density and 
temperature is readily prescribed in terms of R, the distance to the disk axis 
and z, the height above the disk midplane for a given accretion rate. Ma and 
mass of the central star, M*. The density of the disk, n, is given as: 

n(i?, z) = n{R, 0) exp{-z/H{R)) (21) 

with the density in the disk mid-plane being, 

niR, 0) ^ 9.02 X 10^ cm- ( -^^] (^) (^) 
^ ' 1^10-8Moyr-iy \Mq ) M] ) 

(22) 

for R » i?*, as it is the case, and assuming that the main component of 
the disk mass is molecular hydrogen. The scale height of the disk, H{R), is 
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given by, 

H(R) = 3.1 X 10^° cm I 1 ( — ) ( ^\ ] (23) 

^ ' ' lO-8M0yr-M V^0/ VO-1 AU/ ^ ^ 



The disk temperature at a given radius, R, is derived from the a-disk pre- 
scription to be, 

" 0.48f,,^)"'Vi|.V'7^V" (24) 



lO^K llO-SMgyr-i/ \MqJ ^0-1 AU 



5.2 UV irradiation of a-disks 

5*19 varies several orders of magnitude within the dust disk and depends on 
the distance to the star and the opacity of dust grains to the UV radiation. 
UV radiation penetrates through the disk atmosphere creating a photo disso- 
ciation region that has been extensively studied (see i.e. |I2]). It is expected 
that the ionization front driven by the UV radiation splits the atmosphere 
into two regions, an upper ionized layer where H I is fully dissociated and a 
transitional region where the ionization fraction of H I decreases to negligi- 
ble values and the electron density is controlled by Mg I ionization. Further 
inside the disk only X-ray photons are able to penetrate and ionize the dust 
grains. 

Accretion disks around TTSs do not reach the stellar surface; the stellar 
magnetic field disrupts the inner area of the disk. Typically accretion disks 
are considered to extent from some 3-4 stellar radii to some 500 AU (see i.e. 
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however, the inner most area of the disk is devoided of dust grains. The 
UV radiation from the star melts the dust grains to a radius of ~ 0.1—0.2 AU. 

radius, 5*19 = 0.073 and decreases steadily to 6*19 = 2.9 x 10~^ at 500 AU in 
the outskirts of the disk. 

The stellar flux is absorbed by the ionized recombining upper layer and 
by dust absorption. The transfer of radiation the stellar UV flux, in 
this layer is determined by, 



"^^""l^'^^ = -«Bn(i?, zf - n{R, z)auvFuv{R, z) (25) 



where is the case B hydrogen recombination coefficient that at lO'^K is 
as — 2.53 X 10^"*^^ cm^ s~^ [2^, auv is the dust grains absorption cross 
section per hydrogen nucleus, typically auv = 0.61 x 10~^^ cm~^ [2] and 
n{R, z) is the hydrogen density at cylindrical coordinates (i?, z) in the disk. 
The general solution to this equation is: 



Fuv{R, z) = exp —t{z){Feuv,o — / dz'aBn^{R, z) exp —t{z') 

Jo 



(26) 



with r{z) = auv Jq n{R, z')dz' = auvNu being Nh, the column of H I and 
Fuv,o the radiation field at the disk surface. However, the stellar radiation 



^At thermal equilibrium, a dust grain irradiated by a TTS with luminosity L* is heated 
to a temperature T, such that 7rr-g^^^„ ^^^2 = ^'^Irain^'^^ ■ -^^"^ ^ about 2/3 of the sihcates 
melting temperature (T = 2005K), the value of R derived is 0.1 AU, for a solar luminosity 
star. 
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ought to cross the wind departing from the disk that may effectively shield 
the disk against photoionization at early stages (see [6] for an analysis of 
the effect caused by the fast centrifugally driven magnetic winds and [T2] for 
an analysis of the effect caused by slow, thermally driven winds). In such a 
case, the dominant absorbing component is the recombining H I; an optical 
depth of 1 is reached for an H I column of 10^'' — 10^^ cm~^. Moreover, 
the properties of dust grains in the disks atmospheres are ill determined but 
significant deviations from the standard dust grains size distribution in the 
interstellar medium are expectecj^; large dust grains are expected to be more 
numerous due to the high densities in the disk. 

In summary, in the absence of well determined properties of the absorp- 
tion of the UV flux on the disk atmosphere, Sig can be parametrized as, 
5*19^0 exp""^, where S'lg^o is the flux at the disk surface and exp""^ parametrizes 
the extinction as a function of an unknown optical depth, r. Notice that the 
gas and dust columns are very high and that UV radiation will affect only the 
upper atmospheric layer. Typically, the dust grains absorption cross section 
per hydrogen nucleus, ar/v, is ~ 0.61 x 10^^^ cm~^ [2]. Thus, in the innermost 
region of an a-disk, the dust column reaches (l/auv) already at 2; ~ 5H (for 
a standard accretion rate of 10~^MQyr~^). 

^Following 17 , the dust particles size distribution in protostellar disks can be modelled, 
A^d(rgrain) = Nd^r'^L for Vai„ G [0.005 - 

Nd{rgre,in) = ^d^O^tin ' fo^' Vain G (1 " 5]/ini 

According to this distribution, the average dust grain size is < i? >= 0.008/im and the 
maximum charge that such grains can achieve is < Qmax >— 220e^. 
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5.3 Relevance of photoelectric charging in the atmo- 
sphere of an a-disk 

Some hints on the relevance of photoelectric charging of dust grains in the 
accretion disks atmospheres can be derived from a simple evaluation of 
and $g at 2; = 5if for our fiducial dust grain of radius 10 nm. 

On the one hand, the collisional charging term, can be derived using 
as fiducial density n{R, 5H) and deriving the thermal velocity of the gas from 
the temperature of the accretion disk T(R) in the a prescription. Henceforth, 

electrons / T \ ^^V ^ /lO^ K\\ , , 

where Xe is the fraction of electrons in the atmosphere. 

On the other hand, the photoelectric flux, $g, is given in terms of the 
stellar UV radiation flux entering the disk atmosphere, and the optical 
depth, r, as, 

2x 10^0 ^i^^^5i9 exp-- ( ' iQ-^-^^^Q (28) 
cm^'s ' VIO nm/ 

A simple illustration of the competition of both contributions, collisional 
charging versus photoelectric charge loses is shown in Figj6] for a 10 nm 
radius dust grains. In the top panel, $g has been evaluated at the inner 
radius of the dusty disk, 0.1 AU, and plotted as a function of the grain 
charge. The rapid exponential decay is clearly seen. In the bottom panel. 
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the stellar radiation has been assumed to be reduced by a factor 1/e, as would 
correspond to z — 5H. On top of that, the function of $g((5,T) has been 
plotted assuming in both cases that the ionization fraction is Xe — 0.01, and 
that the sticking coefficient is — 1- been plotted for T = 500, 1000, 

1200, 1500 and 5000 K. The final charge of the dust grain is determined by 
the intersection between and the ^'^{Q,T) corresponding to the plasma 
temperature. As the temperature rises, the slope of ^"(Q, T) rises, and the 
equilibrium becomes at a larger charge for the same number of photoionizing 
photons, 5*19. If 5*19 decreases, the region of the exponential $g curve sampled 
by the same family of ^"(Q, T) curves, is at smaller values of Q. 

Notice that the temperature of the gas acts as a regulator of final dust 
grain charge and determines the region dominated by the photoelectric charg- 
ing, as much as the photoionizing radiation flux. 

6 Conclusions 

The effect of the various UV radiation fields on dust grain charging has been 
analyzed by means of numerical simulations to show that TTSs spectra are 
more effective in dust grain charging than the often applied UV-background. 
The photoelectric yield has two main populations, a low energy population of 
electrons with E = 2.5 eV and a high energy population with E = 5.5 — 6 eV; 
the energy dispersion within the populations is ~ 1.3 eV (T ~ 1.5 x 10^ K). 
The high energy population is susceptible of dissociating the H2 and ionizing 
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some low ionization potential species, such as the Mg. These results add 
an additional role to dust on the chemistry of the layers just below the H2 
photoionization front. 

Finally a simple expression has been derived for the photoelectric cur- 
rents from dust grains in TTSs environments that can be introduced on the 
modeling of protostellar disks and a simple example of the relevance of pho- 
toelectric charging is presented for a rather naive a-disk model. 
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8 Appendix 

Velocity Verlet, to start the simulation the initial position, r(to), and velocity, 
v(i(:o), are required. Following must the acceleration, a(t), be determined for 
every time step to progress the system. 



r(t„+i) = r(t„) + v(t„)At + a(t„)At2 



(29) 



(30) 



(see appendix in [25] for more details) 
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Figure 1: Yield as a function of photon energy. Grey curve is model by [29] 
for the decrease at high photon energies, bold black curve is the new proposed 
approximate model, Equ. ID 
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Figure 2: Distribution functions of the photon counts and the cumulative 
value for ionizing radiation. Solid line is the spectrum for a T Tauri star and 
the dashed line is the BB-50kK that both have been normalized to enclose 
an unit area. Grey shaded areas are the cumulative percentage of the spectra 
below the specific wavelength. 
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Figure 3: Charging profiles for a 30 nm grain wlien exposed to a spectrum 
from a T Tauri star (TTS), a Blackbody spectrum at a temperature of either 
10,000 K (BB-lOkK) or 50,000 K (BB-50kK), and the UV background (UV- 
BG). It should be noticed that for the BB spectra linear profiles are obtained 
whereas the TTS results in a S-shaped profile and the UV-BG saturates at 
a plateau level. Inset, charging profiles normalized with the grain radius for 
different grain sizes when exposed to the TTS spectrum. 
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Figure 4: Charging model for an arbitrary sized grain. The fourth-order 
polynomial used to model the dependence between time and the charge state 
of a grain shown in a double logarithmic plot togetether with simulated data, 
the parameters used dbTQ clS stated in Equ. [TTl It should be noted that the 
asymptotic plateau value used in the model is set to 275 e. Insets show the 
simulated data and the fitted model for the charge state as a function of 
time. 
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Figure 5: Energy distribution of the ejected photoelectrons. As the charge 

of the grain increases, only the high energy electrons are able to reach the 
boundary of the system and escape the grain-core. 
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Figure 6: The photoelectric flux leaving a 10 nm dust grain irradiated by 
a TTS spectrum is shown as a function of the grain charge (thick line). 
This current produces a positive charging of the grains to a limit of 275 
electrons. Electrons from the environment are attracted by the charged dust 
grain and an impinging electric current is created. This collisional current 
depends on the electric potential of ^4e grain and on the thermal velocity 
of the electrons. The collisional charging function is represented for various 
thermal velocities of the electrons corresponding to temperatures: 500, 1000, 
1500, 5000 K with decreasing slopes. In the left panel, the photoelectric flux, 
<l>g is represented for 6*19 o = 0.073. In the the right panel, the photoelectric 
flux, $g is represented for Sig^/e = 0.026 



